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ABSTRACT 


A  rapid  pneumatic  transfer  device  has  been  designed  and  constructed 
for  use  with  the  U.S .  Naval  Postgraduate  School  AGN  201  Reactor,  This 
device  permits  rapid  removal  of  the  sample  from  the  Renter  of  the  core  of 
the ^reactor  to  a  shielded  counting  position  to  permit  measurement  of  the 

half “lives  of  short-lived  nuclides.  Preliminary  measurements  of  the  decay 

c  •  x-  r  „r!78m  i  TT£179m  ,  _l£.178m  lir179m 

of  a  mixture  of  Hf  and  Hf  were  made.  The  Hf  and  Hf  nuclides 

were  produced  by  the  neutron  activation  of  a  Hafnium  sample.  Preliminary 

measurements  of  the  resulting  decay  were  made  by  detection  of  the  gamma 

emission.  Decay  data  were  obtained  with  the  use  of  a  scintillation  counter 

and  a  multi-channel  analyzer.  The  measured  half-life  of  Hf 1  ^  was  18.2 

seconds  compared  to  the  published  value  of  19  seconds.  The  half -life  of 
1 7  8m 

Hf  has  not  been  determined  because  of  limitations  of  the  analyzer. 

Proposed  methods  to  improve  measurements  of  both  of  these  half-lives  are 
discussed.  Computer  programs  were  devised  for  least  squares  fitting  of  the 
decay  of  an  independent  mixture  of  two  nuclides  and  for  single  nuclide  decay. 
These  programs  are  described,  and  the  FORTRAN  statement  for  the  decay  of  a 
mixture  of  two  nuclides  is  included. 
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for  his  assistance  and  encouragement  during  the  course  of  this  undertaking. 
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King,  H.E.  Handler,  and  D.G.  Williams^  reactor  operator  F.M.  Spreckhardt, 
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1.  OBJECTIVE 


The  objectives  of  this  thesis  were  the  following 

lc  To  design  and  construct  a  rapid  pneumatic  tube  to  be  used  in 

it 

the  AGN-201  reactor  to  permit  the  rapid  removal  of  an  ir¬ 
radiated  sample  from  the  core  region  of  the  reactor  to  a 
shielded  counting  position. 

2o  To  utilize  this  pneumatic  tube  to  measure  half-lives  of 
short-lived  nuclides  activated  by  neutron  capture. 


The  AGN-201  is  an  experimental  thermal  reactor  with  a  core  which  is  a 
homogeneous  mixture  of  20%  enriched  Uranium  Dioxide  and  Polyethelene • 

The  Reactor  is  graphite  reflected  and  has  a  maximum  flux  of  approximately 
5  x  10^  neut/cm^-sec  at  the  maximum  authorized  power  of  1  KoW. 


1 


II*  TRANSFER  SYSTEM 


A.  General 

In  order  to  determine  half-lives  which  are  in  the  order  of  seconds 
a  procedure  must  be  used  that  allows  the  decaying  nuclide  to  be  counted  al¬ 
most  immediately  after  its  removal  from  the  source  of  activation*  For  our 
measurements  this  was  accomplished  by  the  use  of  a  rapid  pneumatic  tube* 

The  tube  contained  a  sample  holder  or  "rabbit"  which  was  alternately  trans- 
fered  from  an  activation  position**  which  was  the  core  center  of  the  reactor 
(maximum  flux  region) 5  to  a  shielded  counting  position  approximately  six 
feet  outside  of  the  reactor  (Figure  1) *  The  movement  of  the  rabbit  from 

activation  to  counting  position  was  accomplished  by  the  release  of  high 
* 

pressure  air  thru  a  quick-opening  solenoid  valve.  The  time  of  travel  was 
approximately  one-tenth  of  a  second.  When  necessary  the  rabbit  was  returned 
to  the  core  region  for  further  activation  o£  the  enclosed  sample  by  evacu¬ 
ating  the  tube*  This  brought  the  rabbit  back  onto  its1  seat  at  the  acti¬ 
vation  position.  The  cycle  could  be  repeated  many  times. 

Bo  Pneumatic  Tube  and  Accessories 

The  pneumatic  tube  was  constructed  from  7/8"  outside  diameter 
aluminum  alloy  tubing  having  a  wall  thickness  of  *0^9" (Federal  Specification 
WW-T-778b) .  Since  maximum  available  length  of  this  tubing  was  twelve  feet 
and  a  twenty-foot  length  was  needed  to  extend  through  the  reactor,  two  pieces 
were  j  ’red  oj  ’ag  an  internal  sleeve,  which  also  served  as  a 

seat  for  the  rabbito  The  tube  was  inserted  through  the  glory  hole  of  the 
reactor  and  the  rabbit  was  positioned  by  the  seat  at  the  point  of  maximum 
flux  at  the  center  of  the  core. 

Auxiliary  equipment  attached  to  the  aluminum  tube  consisted  of  five 
normally-closed  solenoid  valves  which  controlled  pressure  and  vacuum  in 
the  tube,  and  a  brass  stopping  seat*  The  seat  is  shown  schematically  in 
Figure  2  and  pictured  in  Figure  3.  By  cycling  the  valves,  the  rabbit  was 
moved  from  irradiation  position  to  counting  position  and  back  again* 

Cycling  was  accomplished  manually  with  a  set  of  five  switches  located  on 
the  main  instrument  panel  in  the  reactor  room*  This  panel  also  contained 
pressure  and  vacuum  gauges  for  the  system* 
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vac uufii  and 
p* assure  gauges 
(instrument  pa.  ^1) 


to  vacuum  26  in.  Hg 


Valves; 

V"  1 ,  V 2 ,  V-3,  V~4  Skinner  miniature 
two  way  solenoid  valv  s,  C  serV 
No.  C2DA10S2 

V~5  Skinner  two  way  high  flow  solenoid 
valve,  L  series  No.  L2DB513Q 


PNEUMATIC  TUBE  (Top  View) 


Figure  1.  Pneumatic  Tube 


SCALE:  FULL 

All  dimensions  arc  in  i it r- 1  e 


FAST  PNEUMATIC  TU 
ONE  WAY  VALVE  and  SEAT  ASS;-  BLY 


Figure  2. One  Way  Valve  and  Seat  Assembly 


With  reference  to  Figure  1  a  complete  cycle  of  the  rabbit  is  described. 

The  rabbit  is  placed  into  the  tube  manually  at  the  counting  position  by 
temporarily  removing  the  brass  seat  assembly  at  the  forward  end  of  the  tube. 
Valve  V-4  is  opened .  This  places  a  vacuum  throughout  the  entire  rear  por¬ 
tion  of  the  tube  since  valve  V~3  is  closed  The  rabbit  will  move  somewhat 
but  will  not  enter  the  reactor.  This  is  because  of  an  equilization  of  pres¬ 
sure  in  the  forward  portion  of  the  tube  after  the  rabbit  starts  moving  to- 
wards  the  core*  To  place  the  rabbit  in  the  irradiation  position  (core)  at 
a  pre-set  time,  valve  V~1  is  opened  to  the  atmosphere „  This  creates  a  pres¬ 
sure  differential  of  approximately  one  atmosphere  in  the  forward  portion  of 
the  tube  and  sends  the  rabbit  back  onto  its9  seat0  Valve  V-l  is  then  closed 
and  valve  V-3  is  opened  in  order  to  place  a  vacuum  in  the  entire  tube*  After 
irradiation  of  the  sample  (which  is  contained  in  the  rabbit)  it  is  now  de¬ 
sired  to  move  the  rabbit  to  the  counting  position.  First  valves  V~3  and 
V-4  are  closed  and  then  valve  V~5  is  opened,  releasing  an  air  pressure  of 
thirty  pounds  per  square  inch  against  the  rear  of  the  rabbit.  This  sends 
the  rabbit  very  rapidly  forward  and  into  the  counting  position.  With  the 
rabbit  in  this  position,  valve  V-2  is  opened  for  a  short  time  to  reduce  the 
pressure  in  the  tube  to  atmospheric.  The  cycle  is  then  ready  to  be  repeated 
when  necessary. 

An  accessory  piece  of  equipment  used  in  this  design  was  a  1.5  cubic 
foot  welded  steel  air  pressure  tank  tested  to  120  psi  (Figure  4),  This  was 
connected  to  an  existing  air  compressor  housed  outside  of  the  reactor  build¬ 
ing.  A  pressure  reducing  valve  and  an  air  filter  were  installed  in  the  1/2 
inch  line  between  the  compressor  and  this  tank.  One  end  of  valve  V-5  (quick 
opening)  was  connected  directly  onto  the  tank.  The  opposite  end  was  fitted 
over  the  rear  end  of  the  aluminum  tube  by  means  of  an  *0!  ring  connection 
held  with  three  set  screws.  The  vacuum  system  consisted  of  a  Centro-Pressovac 
4  Pump  driven  by  a  1/4  horsepower  motor  connecting  a  vacuum  reservoir  of  314 
cubic  inch  volume.  All  lines  servicing  pressure  and  vacuum  to  the  aluminum 
tube  through  the  solenoid  valves  were  1/4  inch  outside-diameter,  double 
strength  copper  tubing.  Standard  brass  fittings  were  used  at  all  connections. 

Electrical  connections  between  the  solenoid  valves  and  the  individual 
switches  on  the  instrument  panel  were  made  by  means  of  115-volt,  two-strand, 
copper  wire. 


6 


7 


*  AJfj K  Si<U 


There  were  two  t  yp  s  <  searing  a  r  r  ang-men  t  s  for  the  rabbit  m  the  tube. 
As  previously  mentioned,  the  rabbit  was  stopped  ia  the  irradiation  position 
by  means  of  a  sleeve  inserted  in  the  aluminum  t  ibrig  at  the  desired  place. 
This  also  served  to  strengthen  the  welded  joint.  However  the  seat  (or 
stop)  that  was  located  in  the  counting  position  figure  2)  t  contained  an 
adjustable  spring-loaded  valve  which  permitted  only  one-way  flow  of  air  from 
the  tube  to  the  outside.  The  operating  condition  of  this  valve  was  set  by 
adjusting  the  tension  of  the  enclosed  spring  with  a  threaded  holding  nut. 

This  spring  in  turn  held  a  rubber  gasket  over  the  discharge  opening  of  the 
valve.  The  valve  was  capable  of  passing  air  through  the  opening  at  a  dis¬ 
crete  rate  determined  by  the  tightness  of  the  holding  nut  and  by  the  pres¬ 
sure  differential  between  the  air  in  the  tube  and  the  atmosphere.  The 
necessity  for  this  type  of  valve  was  the  slight  build-up  of  air  pressure  in 
front  of  the  rabbity  (since  there  was  not  a  complete  vacuum  in  the  tube  prior 
to  firing).  This  build  up  can  be  gated  through  the  seat  valve  at  such  a 
rate  as  to  allow  some  of  the  air  to  slow  the  rabbit  down  and  prevent  it  from 
striking  the  seat  with  excessive  force.  This  slowing  down  takes  place  in 
the  last  foot  of  travel  and  can  be  adjusted  as  necessary  to  give  variable 
slow  down  times. 

Co  Rabbit 

The  rabbit  was  constructed  from  balsa  wood  as  shown  in  Figure  5. 

It  weighed  1.07  grams  including  the  hafnium  sample.  The  latter  was  in 
powdered  form  and  was  contained  in  a  gelatin  medicinal  capsule  within  the 
1/8  inch  diameter  hole  located  along  the  long  axis  of  the  rabbit.  Securing 
of  the  rabbit  was  accomplished  by  a  1/8  inch  plug  of  balsa  wood  inserted  and 
glued  in  the  outer  portion  of  the  hole.  Various  other  holding  positions  for 
the  sample  can  be  devised  dependent  upon  the  physical  properties  of  the  sam¬ 
ple  to  be  activated o, 
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Figure  5.  Rabbi c 
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III.  COUNTING  APPARATUS 


After  the  neutron  activated  sample  was  removed  from  the  reactor  core 
and  seated  in  the  counting  position,  the  counting  apparatus  was  used  to 
make  disintegration  rate  measurements  which  were  in  turn  used  to  determine 
the  half-lives  of  the  activated  sample.  A  scintillation  detector  was  chosen 
because  of  the  high  counting  rates  possible.  The  four-inch  lead  shield 
shown  in  Figure  6  was  constructed  around  the  detector  and  rabbit  seat  to 
eliminate  as  much  of  the  background  as  possible .  Automatic  counting  rate 
measurements  were  made  with  the  128-channel  analyzer  by  recording  counts 
observed  during  successive  equal  preselected  intervals  of  time.  The  alter¬ 
nate  system  measured  the  time  required  for  a  preselected  count  using  a  time 
recorder.  Both  systems  have  provisions  for  a  pulse  height  "window11  setting 
to  record  only  pulses  corresponding  to  a  preselected  gamma  energy. 

The  primary  counting  system  is  shown  schematically  in  Figure  7.  The 

detector  consisted  of  a  two-by- two-inch  Nal  (Tl)  scintillation  crystal 

mounted  on  a  DuMont  6292  Photomultiplier  tube.  This  assembly  is  certified 

by  the  manufacturer  to  be  capable  of  7.6%  resolution  (defined  as  the  energy 

width  at  half  maximum  of  the  photopeak  divided  by  mean  photopeak  energy)  of 
137 

the  .662  MEV  Cs  gamma.  The  scintillation  counter  was  used  with  an  ex¬ 
ternal  high  voltage  power  supply  to  replace  the  internal  high  voltage  power 
supply  (in  the  RCLiac  128~channel  analyzer)  which  was  found  to  have  consider¬ 
able  drift.  For  measurement  of  counting  rates  with  this  system  the  "gross 
mode"  of  analyzer  operation  was  used.  The  number  of  counts  observed  during 
each  of  125-constant  preselected  time  intervals  were  photographically  recorded 
from  the  visual  display  of  the  analyzer.  The  counts  corresponded  to  a  given 
gamma  energy  range  (determined  by  the  window  setting).  The  proper  window 
setting  was  readily  determined  by  using  the  1 28 ~ channel  "pulse  height  analyzer 
mode"  of  operation.  (See  calibrations  page  16). 

An  alternate  counting  system  is  shown  schematically  in  Figure  8.  It  in¬ 
cludes  the  scintillation  crystal  described  above 9  a  cathode  follower,  and  a 
non- over loading  linear  amplifier  with  provisions  for  pulse  height  window 
settings.  The  decimal  scaler  has  been  modified  to  deliver  pulses  at  100  or 
1000  count  intervals  to  a  two  channel  recorder.  The  second  channel  of  the 
recorder  was  used  for  recording  time  calibration  pulses.  The  high  voltage 
for  photomultiplier  operation  was  supplied  by  a  regulated  high  voltage  power 
suppy.  The  occurrence  of  each  100  or  1000  counts  was  recorded  as  a 
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Figure  6*  -ead 


COUNTING  SYSTEM  USING  128  CHANNEL  ANALYZER 


Figure  7.  Counting  System  Using  128~Channel  Analyzer 


ALTERNATE  COUNTING  SYSTEM 
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Figure  8.  Alternate  Counting  System 


tick  mark  on  the  recorder  chart  which  was  moving  at  a  constant  rate.  The 
elapsed  time  between  two  marks  was  determined  by  measuring  the  separation  on 
the  chart  paper.  This  time  interval  was  compared  with  the  time  calibration 
made  on  the  other  channel.  The  equipment  was  assembled  and  tested,  but  has 
not  actually  been  used  to  obtain  data. 

Figure  9  is  a  photograph  of  all  of  the  assembled  apparatus. 
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Three  calibrations  were  required  for  the  measurement  of  half-lives  using 
the  128  channel  analyzer;  the  live  time  and  dead  time  calibration  for  the 
rrgross  mode1*  of  operation*  and  the  energy  (pulse  height)  vs.  channel  number 
calibration  for  the  npulse  height  analyzer11  mode  of  operation, 

The  live  time  is  the  time  available  for  counting  pulses  in  each  of  the 
channel  So  The  live  time  is  nearly  equal  to  the  time  per  channel,  differing 
only  by  the  constant  time  required  for  the  analyzer  to  shift  between  channels. 
The  time  per  channel*  referred  to  as  At,  cannot  be  readily  measured  directly, 
and  thus  it  is  necessary  to  refer  to  the  live  time,  which  closely  approxi¬ 
mates  At,  for  an  estimate  of  the  statistical  variation  of  At,  The  cali 
bration  was  accomplished  by  inserting  pulses  from  a  10  KC  crystal  controlled 
oscillator  into  the  second  stage  of  the  photomultiplier  preamplifier.  An 
analysis  of  the  time  calibrations  made  during  an  eight  day  period  gave  live 
times  of  „1902  +  o0006  seconds,  O3804  j-  .0007  seconds,  and  7601  +  0005 

seconds  for  the  .19  x  1,  „19  x  2,  and  ,19  x  4  seconds^per  channel  time 

settings.  The  ^t's  for  the  same  settings  are  <,192*  0384,  and  0  768  seconds 
per  channel.  The  live  time  drift  observed  during  the  eight  day  period  ap¬ 
peared  to  be  within  the  above  statistical  variations  of  the  live  times.  As 
an  additional  safeguard,  checks  of  live  times  were  made  both  before  and  after 
the  taking  of  data. 

Attempts  to  determine  the  dead  time  -  the  time  following  the  incidence 
of  a  pulse  that  the  counter  will  not  detect  a  second  pulse  -  were  unsuccess¬ 
ful.  Measurements  were  taken  using  the  two  source  method,  which  gave  dead 
times  of  150,  190,  and  170  microseconds  for  three  successive  trials  for  three 
different  energy  settings.  This  variation  with  energy  is  not  anticipated 
from  the  design  of  the  analyzer „  The  dead  time  measurements  were  made  some 
time  after  the  decay  measurements*  and  it  is  suspected  that  instrument  mal~ 
function  at  the  time  of  calibration  is  responsible  for  the  poor  results.  If 
the  actual  instrument  dead  time  i s  as  large  as  150/^  sec*  the  maximum  number 
of  counts  that  can  be  recorded  in  a  .19  second  interval  without  excessive 
coincidence  correction  (defined  for  our  purposes  as  5%)  is  300,  This  sets 
a  limit  of  6%  statistical  error  on  the  counting  rate  determination  by  a  sin¬ 
gle  channel.  During  the  calibration  it  was  not  possible  to  record  sufficient 
counts  using  the  ,19  x,  2  second  channel  time  setting  to  permit  statistically 


s out  ^  measurement s  It  was  noted  that  at  appr o* ima tel y  800  counts  per 
channel  the  subsequent  positioning  ot  the  calibration  sources  closer  to 
the  scintillation  crystal  initially  had  no  effect,  Then ^  even  closer 
positioning  caused  the  count  to  drop  to  zero,  This  apparent  saturation 
effect  has  not  been  reconciled  with  the  counting  of  10,000  counts  per 
second  during  the  live  time  calibration  Since  it  is  anticipated  that 
the  alternate  counting  system  will  be  used  in  the  future,  (see  Recommen¬ 
dations)  no  further  attempts  were  made  to  determine  the  precise  dead  time* 
The  energy  versus  pulse  height  calibration  was  made  using  the  "pulse 

height  analyzer"  mode  of  operation  of  the  128-channel  analyzer „  Samples 
137  198  153 

of  Cs  ,  Au  ,  and  Sm  having  known  gamma  energies  of  6616  MEV  for 
Cs,  *4118  and  „069  MEV  for  Au,  and  ,1032  MEV  for  Sm  were  used  to  obtain 
the  calibration  curve  (Figure  10)  a  The.  curve  is  a  plot  of  the  channel 
number  at  which  the  maximum  count  occurs  in  the  photopeak,  versus  the 
known  energy  of  the  photopeak* 
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ANALYZER  CHANNEL  NUMBER 
figure  lO.Gajmna  Energy  Calibration 

18 


V.  RESULTS 

Measurements  of  the  half  lives  of  two  neutron  activated  nuclides  of 

Hafnium  were  attempted  using  the  128-channel  analyzer  and  pneumatic  system 

already  described,  A  computer  program,  described  in  Appendix  1,  was  used 

to  make  a  least  squares  fit  of  the  experimental  data  for  the  19  second  half- 

life  Hf^^™  nuclide  to  a  single  exponential  decay  curve.  Attempts  to  verify 

1 7  8m 

the  4,8  second  half-life  of  the  Hf  nuclide  were  unsuccessful , 

The  Hafnium  isotopes  Hf^  ~  and  Hf^^  have  thermal  neutron  cross  sections 

2 

of  370  and  80  barns  respectively.  When  a  sample  containing  these  nuclides 


is  exposed  to  a  neutron  flux  the  following  reactions  take  place  Hf 


177 


Tir178m  xa  ,  ,Tir178  1 

Hf  -f  q  and  Hf  +  n 

o 


Hf 


179m 


f  o  The  resulting  nuclides 


Hfl/8m  an(j  ^179m  are  metastable  and  decay  by  gamma  emission  to  the  ground 
state  with  half  lives  of  4,8  seconds  and  19  seconds  respectively.  The  decay 
schemes  are  shown  diagramatically  in  Figure  11,  Gamma  energies  emitted  are 
.089,  .093,  .214,  .326,  and  ,427  MEV  for  the  Hf^^m  nuclide  and  .160  and 


,215  MEV  for  the  Hf 


179m  1,3 

nuclide. 


The  sample  used  consisted  of  50  mg  of  Hafnium  in  oxide  form  enriched  to 

59,08%  Hf^7  and  also  contained  31,737o  Hf^  ^  (Appendix  III),  The  sample  was 

8  9  2 

exposed  to  neutron  fluxes  of  10  and  2,5  x  10  neutrons/cm  —second  for  30 
seconds,  (These  correspond  to  reactor  powers  of  2  and  50  watts  respectively,) 
The  128-channel  analyzer  was  initially  used  in  the  pulse  height  analyzer 
mode  of  operation  to  examine  the  gamma  spectrum  of  the  activated  Hafnium, 

Since  both  Hf^^  and  Hf^^  are  present  in  large  percentages,  gamma  energy 
peaks  were  expected  at  ,160  and  ,215  MEV  from  the  decay  of  Hf ^  ^  and  ,089, 
,093,  ,214,  ,326,  and  ,427  MEV  from  the  decay  of  Hf^^m  (Figure  11),  Major 
peaks  were  observed  at  .07  +  ,01  MEV  and  ,22  ,01  MEV.  Smaller  peaks  were 

found  at  .33  +  ,01  and  .44  +  .01  MEV.  The  peaks  were  attributed  to  gamma 
activity  because  of  the  beta  shielding  inherent  in  the  detector  arrangement. 
Observation  of  the  sample  immediately  after  activation,  and  then  after  30 
seconds  decay,  showed  the  ,07  and  .22  MEV  peaks  to  be  due  primarily  to  the 
longer  lived  nuclide  i.e.,  Hf^^m,  The  ,22  MEV  activity  is  attributed  to 
the  19  second  nuclide,  and  the  .07  MEV  corresponds  to  the  Hafnium  X-ray. 

A  window  of  approximately  .04  MEV  width  was  established  at  . 22  +  .01 
MEV  and  the  gross  mode  of  analyzer  operation  was  used  to  obtain  the  decay 
curve  shown  in  Figure  12,  A  At  of  ,19  x  4  seconds  per  channel  was  chosen 
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MEV 


MEV 


DISINTEGRATION  ENERGY  AND  SCHEME 

uf,178m  179m 

Hf  and  Hf 


Figure  11.  Di sintegration  Energy  and  Schemes 


Hf178" 


and  Hf 


179m 
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Figure  12.  Decay  of  Neutron  Activated  Hf  Sample 
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to  meet  the  requirement  that  the  decay  time  is  long  compared  to  the  inter¬ 
val  of  measurement o  (Appendix  I.  explains  this  requirement,)  Analysis  of 
tliase  data  using  the  linear  least  squares  fitting  program  described  in 
Appendix  I,  gave  a  half-life  of  18,15  +  .05  seconds.  A  50  microsecond  dead 
time  correction  was  used.  The  background  was  experimentally  determined  to 
be  less  than  1  count  per  second  and  was  neglected.  The  possibility  of 
systematic  errors  in  the  determination  and  other  experimental  difficulties 
described  below  preclude  the  reporting  of  this  value,  although  the  statistical 
error  would  indicate  an  improvement  over  the  published  value.  The  use  of 
the  alternate  counting  system  should  permit  a  more  satisfactory  half-life 
determination . 

178m 

The  4.8  second  Hf  nuclide  has  reported  gamma  energies  of  .089, 

1  3 

.093,  .214,  .326,  and  .427  MEV.  *  Attempts  were  made  to  isolate  the  .427 
MEV  gamma  using  a  window  setting  of  .43  +  .01  MEV  with  a  .04  MEV  width. 

After  activation  at  70  watts  for  30  seconds  using  a  4t  of  ,19  x  2  seconds 
per  channel  the  decay  curve  shown  in  Figure  13  was  obtained.  The  initial 
portion  of  the  curve  has  a  slope  corresponding  to  the  4.8  second  half-life, 
but  the  "hump11  appearing  later  causes  serious  question  as  to  the  validity 
of  this  value.  Attempts  to  eliminate  the  hump  using  lower  reactor  powers 
were  successful  (it  did  not  appear  at  reactor  powers  less  than  20  watts); 
however,  this  also  resulted  in  insufficient  activity  to  identify  the  half- 
life  with  any  degree  of  precision.  It  was  noticed  in  this  process  that  re¬ 
duction  in  power  or  activation  time  caused  the  hump  to  appear  at  an  earlier 
time.  Measurements  taken  at  slightly  different  window  settings  with  similar 
reactor  power  and  activation  time  do  not  show  the  characteristic  hump,  but 
instead  show  a  9-10  second  half-life.  From  this  information  it  may  be  postu¬ 
lated  that  the  hump  is  caused  by  coincident  ,215  MEV  gammas.  They  initially 
occur  at  counting  rates  sufficiently  high  to  cause  coincidence  of  more  than 
two  such  pulses.  The  resulting  composite  pulse  is  rejected  because  it  has 
too  high  an  amplitude.  As  the  number  of  .215  MEV  pulses  decrease,  the  pos¬ 
sibility  of  the  coincidence  of  pairs  of  pulses,  as  opposed  to  the  coincidence 
of  more  than  two  pulses,  increases  and  thus  pulses  equivalent  to  .430  MEV 
gammas  are  counted.  The  number  of  such  pulses  decreases  with  a  half-life 
equal  to  1/2  the  half-life  of  the  decaying  nuclide  i.e.,  9,5  seconds.  It  is 
hoped  that  by  the  selection  of  a  different  gamma  energy  and  use  of  the  alter¬ 
nate  counting  equipment  that  it  will  be  possible  to  resolve  the  short-lived 
activity. 
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VI ,  RECOMMENDATIONS 


Several  modifications  and  additions  to  the  present  pneumatic  tube 
system  are  recommended  in  order  that  the  system  be  of  maximum  benefit  to 
the  LLS.  Naval  Postgraduate  School  in  future  work.  These  are  outlined  as 
follows „ 

1.  Increase  the  vacuum  reservoir  tank  from  the  present  314  cubic 
inch  volume  to  approximately  1500  cubic  inches 0  This  will  permit  a  faster 
evacuation  of  air  from  the  tube  prior  to  transfer  of  the  rabbit  into  the 
core  position.  This  allows  a  second  rapid  cycling  of  the  rabbit  if  desired. 

2.  Install  an  automatic  variable  timing  mechansim  to  replace  the 
present  manual  switches.  (The  switches  effectively  determine  pressure  and 
vacuum  in  the  tube  by  controlling  the  five  solenoid  valves.)  This  would 
make  the  cycling  process  a  completely  automatic  procedure  with  the  capability 
of  changing  irradiation  and  counting  times  as  necessary. 

3.  Construct  a  special  seating  arrangement  at  the  counting  position 
of  the  pneumatic  tube.  This  feature  would  permit  the  removal  of  the  rabbit 
from  the  counting  position  without  removing  the  seat  from  the  pneumatic  tube. 
In  constructing  this  modification,  thought  should  be  given  to  the  possibil¬ 
ity  of  also  measuring  beta  activity  at  the  counting  position. 

The  following  recommendations  are  made  concerning  the  measurements  of 
short  half-lives: 

1.  That  the  alternate  system  of  counting  (Figure  8)  be  used  in 
lieu  of  the  128-channel  analyzer.  This  would  result  in  a  reduction  of  both 
systematic  errors  and  statistical  errors.  However,  the  128-channel  analyzer 
is  very  useful  for  exploratory  decay  measurements  and  measuring  the  energy 
spectrum  of  a  particular  sample. 

2.  That  measurements  of  Hf  continue  using  the  alternate  system 
of  counting. 

3.  That  measurements  be  made  on  Yb^^  and  Ba^*\  50  mg  samples 

176  135 

of  Yb  and  Ba  enriched  in  Yb  and  Ba  respectively  are  presently  available. 
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APPENDIX  I 


Theory  and  Computer  Program 

Assume  that  the  activity  of  a  particular  sample  of  a  mixture  of  two 
independent  nuclides  having  different  half-lives  and  a  constant  background 
is  given  by  the  expression. 


A. 


n  =  nioe 


+  n2Qe 


A 


2t 


+  c 


(1) 


where  and  are  the  initial  activities  at  t  =  0,  A^  andA2  are 

the  two  decay  constants  and  c  is  the  constant  background  activity.  It  is 
desired  to  experimentally  determine  the  five  parameters  appearing  in 
equation  (1),  and  in  particular  to  determine  the  decay  constants  A  andA2* 


Experimentally  one  measures  the  total  number  of  counts  occurring 

in  a  known  interval  of  time  At.  It  is  common  practice  to  specify  the 
decay  time  t  to  be  at  the  midpoint  of  the  interval  A  t.  Referring  to 
Eq  (1) 


AN  = 


and  the  experimentally  determined  mean  counting  rate  over  the  interval  A t  is 

f  t  * 

AN  \ 


N  = 


At 


Using  Eq  (1) 


At 


cU 


t  - 


At 


N  -  n-tL 

A  At 


-  X  t  /  X  At 


A,  At 


n 


(2) 


n 


2t> 


-A„t 


\At 


Az^t 


-  A,  At 


-+ 


C 
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we  now  make  the  following  assumptions  about  A  t  and  the  decay  constants, 
A  j>  (Nuclide  1  decays  faster  than  nuclide  2C) 


At  < 


(The  measuring  interval  is  small  compared  to  the  mean  life 
of  nuclide  1,  and  thus  even  smaller  than  the  mean  life  of 
nuclide  20) 


Expanding  the  exponential  terms  in  a  power  series,  and  neglecting  terms 
of  smaller  than  second  order  results  in  the  approximation  to  Eq.  (2)  of. 


N  = 


n  *=> 

\0 


(AM)'] 

2^  - 


n,  e 

?.o 


2.4 


+  C 


(3) 


If  At  is  1/10  of  the  mean  life,  then  the  correction  factor  (second  order 

— 4 

term)  in  the  first  bracket  is  ^4  x  10  and  the  corresponding  correction 
for  the  second  term  is  even  smaller.  Henceforward  it  will  be  assumed  that 
A  t  is  always  chosen  small  enough  so  that  the  above  correction  terms  are 

negligible  compared  to  the  statistical  uncertainty  of  the  experimentally 

•  * 

determined  N.  Thus  N  =  n  and  the  form  of  Eq.  (1)  can  be  used.  If  this 
were  not  so,  Eq.  (2)  or  Eq.  (3)  would  have  to  be  used  for  the  experimental 
determination  of  A^  andA^* 


Therefore  N  =  n^e 


A,  £  .  -  A~t 

L  j-  n  P  £ 


+  n20e 


+  c 


(4) 


where  N  = 


An 

At 


and  t  is  always  specified  to  be  at  the  midpoint  of  the 


constant  time  interval  At. 

If  activity  2  and  the  background  activity  were  completely  missing  then 

N  =  niQe  >  and  ^ n  N  -  Jh\  n^-  A^t  is  the  resulting  linear  equation 

which  can  then  easily  yield  a  least  square  determination  of  A 1 . 
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If  activity  2  is  not  present,,  but  a.  background  activity  c  is  present,  then 


/n  =  ,/n  n^-X^t  can  be  least  squares  fitted  to  determine  A^. 

This  procedure  assumes  that  c  is  previously  experimentally  determined.  If 
c  is  not  independently  determined  then  the  best  values  of  n  ,  and  c 

will  be  determined  by  a  least  squares  fit  to  the  equation 


A, 


N  =  n|Qe  +  c- 


A  least  squares  fit  of  the  observed  counting  rate  N  and  t  has  been 

* 

devised  using  a  Taylor  series  approximation  .  If 

r  -  -AA  .  ~A-d  t 

N  ~  4  ^  is  differentiated  with  respect  to 


^  *  Y  '  "2.0 


d  A, 


nio  / 

—  c 
t 


then 

A,  t 


and 

C, 

SN 

-AA 

_ 

-  £ 

e 

Hro  ~  A 

'A* —  '2. 

..  t 

.  and 

c>N 

\ 

1 

A, 

for  a 

particular 

point  n . , 

t . 

A  A 


using 


only  the  first  derivatives  becomes 


nt  ^  n  V. 

l 


AtL 


-  "is  -  VXAX. 


t 


n 


--  £ 


^  AX,  ^  Ac  (5) 


where  A,  ,\2,n,  .  A  A  ^  c  represent  approximations  to  \  )A2o,  n^,C 


and  ZAn,  ^  r>  A 

•O  »o 


lo 


,  Anld  ■  AX,  X,  A,  ,  A  A, *  X,- X, }  Ac  =  c  - 


c  - 


The  best  values  of  ;nin,  n  ,  A,,  A,  and  c  are  found  first  by  writing  i 
*  1U  2U  1  A 


This  method  has  been  suggested  by  Prof.  G .W .  Rodeback  who  has  applied  it  to 
least  squares  fitting  of  hyperbolic  sine  and  other  non-linear  functions. 
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equations  of  the  form  : 


AX. 


(6) 


Then  the  corrections  A  h  ,  aX  aA  ,  and  Ac  are  found  by  matrix 

methods  resulting  in  improved  values  forA-^^A^j  h^Q,  ^203  C#  This  Proce 
is  repeated  until  the  corrections  become  arbitrarily  small.  In  each  step 
the  term 

-  H 


£l  :(N.  A  J 


is  minimized,  the  weighing  factor,  is  inversely  proportional  to  the 

square  root  of  the  number  of  counts  corresponding  to  the.ith  interval  since 

o 

this  represents  the  standard  deviation  of  NL  . 

The  program  for  the  IBM  1604  Computer ,  as  devised  by  personnel  of  the 
U.S.  Naval  Postgraduate  School  computer  center,  is  attached.  Results  from 
trial  runs  using  this  method  with  both  numerically  devised  data,  and  experi¬ 
mentally  determined  data  for  silver  activated  in  the  AGN  201  reactor,  have 
been  obtained.  In  the  case  of  silver  activities,  the  half  lives  were  de¬ 
termined  as  23.2  seconds  and  2.3  minutes,  which  is  in  agreement  with  the 
reported"*  values  of  24  seconds  and  2.3  minutes  for  Ag '  ^  and  Ag^^  re¬ 
spectively. 
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JOB  WILLIAMS 

PROGRAM  FLUX 

ODIMFNSION  X! 150,5) ,PHIO(  150) , PHIC!  I  SO  )  ,  T  (  1  SO )  ,  DELPHI  I  (  1  50  )  ,  W  (  1  SO  ) , 
1H(5),XW(5, 150), A (5,6) ,DELC0(5) 

33  READ  1,  C1.CK1.C2.CK2.C3 

1  FORMAT! 5F 10.0 ) 

KK  =  0 

READ  2,N,DELT,L 

2  FORMAT! 13. F3. 1,13) 

READ  31  , ! T! I ) ,PHIO( I ) , I- 1 ,N) 

31  F0RMAT(6(F4.2,3X,F5.0)  ) 

DO  32  I  =2 , N 

32  T!I )  =  T ! I )  +  T( I-  1  ) 

DO  200  1=1, N 

CHANL  =  L+I-l 
CNST=38.E-6  ♦  0 . 5E-6  *CHANL 
200  PHIOt I )  =  PHIO(  I )/!  1 .-PH  10! I ) *CNST /DEL  T ) 

DO  100  I  =  I  ,  N 
100  W! I )  =  SORTF! PHIO! I  )  ) 

20  DO  10  1=1, N 
T 1  =  T (  I  ) 

EXIT  =  EXPF(CKWTI) 

EK2T  =  EXPF ! CK2*T I  ) 

PHIC! I )=C 1*EK 1T+C2*EK2T  +C3 

X! I ,  1  )  =  EK IT 

X! I, 3)  =  EK2T 

X!  1,2)  =  C1*TI*EK1T 

X!  1 ,4 )  =  C2  *T  I  *  EK.2T 

X!  I , 5 )=  1  .0 

10  DELPHI (  I )  =  PH  10(1)  -  PHIC!  I  ) 

PRINT  8  (DELPHI!  I  )  ,  I  -=  I.N) 

DO  11  J=1,  N 

W=W(  J  ) 

DO  11  1=  1,5 

11  XW ( I , J  )  =  W*X ( J , I  ) 

DO  13  1=1,5 

H( I )=0. 

0013  K=  1  ,  N 

13  H(  I  )  =  H{ I )  +XW! I ,K)*DELPHI(K) 

DO  15  I  =  1,5 

DO  15  J  =  1,5 
P=0. 

DO  14  K=  1 , N 

14  P=P  +  XW{ I ,K)*X(K,  J  ) 

15  A! I, J)=P 

DO  12  1=1,5 

12  A! I ,6)  =  H( I  ) 

PRINT  300 

300  FORMAT (4H  DEL) 

CALL  J0RDAN2(A,5,DELC0) 

DNORM  =  0. 

DO  16  I =  1 ,  N 

16  DNORM  =  DNORM  +  DELPHI ! I )*OELPHI  (I) 

E  ST  =  0 . 

DO  17  1=1,5 

17  EST  =  EST  +  DELCO! I ) *H(  I  ) 

PRINT  301  ,  DNORM, EST 

301  FORMAT! 2E 12.5) 

PRINT  7  KK,C I ,CK  1  ,C2,CK2,C3 
7  F0RMAT(2X, I2,5X,5t 15.6) 

PRINT  8  (DELCO! I ), 1=1,5) 


/ 


8  FORMAT (  7X ,  5E17.7) 

Cl  =  C  1  +DELCO l  I ) 

CK 1 =  CK1  +DELCO I 2 ) 

C2  =  C2  +  OELCO ( 3 ) 

CK2  =  CK2  ♦  DELC014) 

C3  =  C3  +  OELCO ( 5 ) 

KK  =  KK  +  1 

I F (  KK  -20)  20,20,30 
30  PRINT  40 
40  F0RMAT(5H1J0B2) 

GO  TO  33 
ENO 

SUBROUTINE  J0RDAN2 ( A , N , X ) 
DIMENSION  A ( 5 , 6 )  i  X(5) 

K  =N+  1 

11  IF  ( K—  1  ) 13,6,  15 
15  0=0. 

00  2  r=2,K 

IF  (ABSF{A{ I —  1  ,  1  ))  -0)  4,4,3 

3  L  =1-1 

2  D  =ABSF( AIL, 1 ) ) 

4  IF ( L- 1 )5, 6, 5 

5  00  7  J=  I , K 
D  =  A  (  L ,  J  ) 

AIL,  J)  =  A{  l, J) 

7  Al  1  ,  J  )  =  0 

6  DO  8  I  =  1 , N 

8  XII)  =  A  I  I  ,  l  > 

IF  IK-1  )  12,  13,  12 

12  DO  10  J=2 » K 

0  =  Al  1,J)/X{  1) 

00  9  1=2, N 

9  A I  I-  l , J-  1  )  =  A I  I , J  )  -XII)*0 

10  AIN, J-  1  )  =0 

K  =  K—  I 
GO  TO  l I 

13  CONTINUE 
RETURN 
END 

END 


i 
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APPENDIX  III 

Analysis  of  Hafnium  Sample 


ANALYSIS  OF  HAFNIUM  SAMPLE 
(as  reported  by  Oak  Ridge  analysis  report) 


ELEMENT 

Hafnium 

ISOTOPIC  ANALYSIS 

ISOTOPE 

177 

ISOTOPE 

ATOMIC 

PERCENT 

PRECISION 

SERIES 

GP 

174 

Trace  Only 

SAMPLE 

827(a) 

176 

0.99 

+  0.07 

177 

59.08 

+  0.52 

178 

31.73 

+  0.49 

179 

4.53 

+  0.17 

180 

3.68 

+  0.11 

The  limits  quoted  above  are  an  expression  of 
the  precision  of  this  measurement  only.  The 
error  is  estimated  at  less  than  1%  from  known 
sources  of  systematic  errors. 


SPECTROSCOPIC  ANALYSIS 
ELEMENT  PERCENT 


Cu 

Fe 


0.1 

0.2 


The  spectrographic  results  reported  herein  are  semi-quantitative 
estimates  and  should  not  be  interpreted  or  construed  to  be  pre¬ 
cise  quantitative  determinations. 
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